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The critical behaviors of the bimodal and Gaussian Ising spin glass (ISG) models in dimension 
four are studied through extensive numerical simulations, and from an analysis of high temperature 
series expansion (HTSE) data of Klein et al. (1991). The simulations include standard finite 
size scaling measurements, thermodynamic limit regime measurements, and analyses which provide 
estimates of critical exponents without any consideration of the critical temperature. The higher 
order HTSE series for the bimodal model provide accurate estimates of the critical temperature 
and critical exponents. These estimates are independent of and fully consistent with the simulation 
values. Gomparisons between ISG models in dimension four show that the critical exponents and the 
critical constants for dimensionless observables depend on the form of the interaction distribution 
of the model. 
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INTRODUCTION 


Renormalization Group Theory (RGT) for thermody¬ 
namic phase transitions [l[ and the Edwards-Anderson 
model for Ising Spin Glasses (ISGs) 0 were introduced 
almost simultaneously forty years ago. Ever since it has 
been tacitly assumed as self-evident that the standard 
RGT universality rules should apply to ISGs. As far as 
we know there is no rigorous theoretical proof that this 
ISG hypothesis holds, though confirmations have been 
reported a number of times based on numerical data 
[Jli- The universality principle states that for all sys¬ 
tems within a universality class the critical exponents are 
strictly identical and do not depend on the microscopic 
parameters of the model. All ISG models in a given di¬ 
mension are supposed to be in the same universality class, 
on the assumption that the form of the interaction dis¬ 
tribution is an irrelevant microscopic parameter. 

Thus in the family of simple ferromagnets, within a 
universality class of models having space dimension d and 
spin dimensionality n, all models have identical critical 
properties corresponding to an isolated fixed point in the 
renormalisation group flow. However, diluted ferromag¬ 
nets of given d, n have a different (dilution independent) 
set of critical exponents, with values which correspond 
to a separate isolated fixed point Bi- For a few spe¬ 
cial cases of spin models in dimension two (discussed for 
instance in Ref. i) the critical behavior is more compli¬ 
cated and corresponds to a line of fixed points rather than 
an isolated fixed point; critical parameters vary continu¬ 
ously according to motion along the line, produced by a 
marginal operator. From the empirical ISG data below 
there appear for the moment to be two possible scenar¬ 
ios : two classes of ISGs (such as models with continu¬ 
ous distributions and those with discrete distributions) 
or alternatively ISG exponents which vary continuously 
with a parameter such as the kurtosis of the interaction 
distribution. In any case it has been stated by authori¬ 


tative authors that “classical tools of RGT analysis are 
not suitable for spin glasses” l3“12| although no explicit 
theoretical predictions have been made so far concerning 
the important question of universality. 

Here we combine numerical simulation and high tem¬ 
perature series expansion (HTSE) data on the bimodal 
and Gaussian ISGs in dimension four so as to obtain ac¬ 
curate and reliable values for the critical parameters in 
this model. We discuss a number of different methods 
for exploiting numerical data, and show that for each 
model these are consistent. Comparisons between these 
and other estimates on ISG models in the same dimension 
but with different interaction distributions show that the 
critical exponents and the critical constants for dimen¬ 
sionless parameters depend on the form of the interaction 
distribution. 


NUMERICAL TECHNIQUES 

The Hamiltonian is as usual 

T~L = ~ JijSiSj ( 1 ) 

ij 


with the near neighbor symmetric distributions normal¬ 
ized to (Jfj) = 1. The Ising spins live on simple hyper- 
cubic lattices with periodic boundary conditions. We 
have studied the bimodal model with a ±J interaction 
distribution and the Gaussian interaction distribution 
model. We will compare with published measurements 
on these and other 4d ISGs. We will use the inverse 
temperature /3 = = I/T, with the normaliza¬ 

tion above, or alternatively w = tanh^(/3), to signify the 
temperature. The spin overlap parameter is defined by 



( 2 ) 
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where A and -B indicate two copies of the same system 
and N is the number of sites. 

The simulations were carried out using the exchange 
Monte Carlo method for equilibration using so called 
multi-spin coding, on 2^^ (up to L = 7) or 2^^ (for larger 
L) individual samples at each size. An exchange was at¬ 
tempted after every sweep with a success rate of at least 
30%. At least 40 temperatures were used forming a ge¬ 
ometric progression reaching down to /3max = 0.55 for 
the bimodal model and /3max = 0.60 for the Gaussian 
model. This ensures that our data span the critical tem¬ 
perature region which is essential for the FSS fits. Near 
the critical temperature the P step length was at most 
0.03. The various systems were deemed to have reached 
equilibrium when the sample average susceptibility for 
the lowest temperature showed no trend between runs. 
For example, for L = 12 this means about 200000 sweep- 
exchange steps. 

After equilibration, at least 200000 measurements were 
made for each sample for all sizes, taking place after 
every sweep-exchange step. Data were registered for 
the energy E{P,L), the correlation length f(/3,L), and 
for the spin overlap moments. In addition the corre¬ 
lations {E{P,L),U{P,L)) between the energy and ob¬ 
servables U{P,L) were also registered so that thermody¬ 
namic derivatives could be evaluated using the relation 
dU{P,L)/dp = {UiP,L)EiPM) - (C/(/3,L))(A(/3,L)) 
where E{P,L) is the energy [1^. Bootstrap analyses of 
the errors in the derivatives as well as in the observables 
U{P,L) themselves were carried out. 


FINITE SIZE SCALING 

ISG simulations are much more demanding numeri¬ 
cally than are those on, say, pure ferromagnet transitions 
with no interaction disorder. The traditional approach to 
criticality in ISGs has been to study the temperature and 
size dependence of dimensionless observables, principally 
the Binder cumulant g(/3,L) and the correlation length 
ratio L)/L, in the near-transition region and to es¬ 

timate the critical temperature and exponents through 
finite size scaling (FSS) relations after taking means over 
large numbers of samples. Finite size corrections to scal¬ 
ing must be allowed for explicitly which can be delicate 
as the range of sizes L is generally small. On this FSS 
approach the estimated values for the critical exponents 
are very sensitive to the critical inverse temperature Pc 
estimates. Here we first obtain estimates for Pc and the 
critical parameters using FSS. 

The data for standard dimensionless observables, the 
Binder cumulant 



FIG. 1: (Color online) The 4d bimodal ISG. FSS for the cor¬ 
relation length ratio ^{P,L)/L at fixed inverse temperatures 
P = 0.5100, 0.5075, 0.5050, 0.5025, 0.5000,0.4975 from top to 
bottom, against 1/L^’^. Red squares : present data, black 
circles : read from Ref. El- Dashed line : critical behavior. 
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FIG. 2: (Color online) The 4d bimodal ISG. FSS for the 
Binder cumulant g{P, L) at fixed inverse temperatures /3 = 
0.5100, 0.5075, 0.5050, 0.5025, 0.5000, 0.4975 from top to 
bottom, against 1/L^’^. Red squares : present data, black 
circles : read from Ref. [T^. Dashed line : critical behavior. 


Figs. [T] and [2] as functions of size at fixed temperatures 
near Pc- We show data from the present simulations to¬ 
gether with data taken from the appropriate figures in 
Ref. E^. The standard FSS expression for a dimension¬ 
less observable U = U{P,L), valid in the critical region 
is 


and the correlation length ratio 5(/3, L)/L are shown in U = Uc + AL + B{P — Pc)L^^’' (l + CL “-!-•••) (4) 
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where Uc = U{Pc,oo). We plot g{P,L) and ^{I3,L)/L 
against l/L^^ with uj = 1.2. The choice of this value, 
close to the estimate w = 1.04(10) of Ref. [l^, will be 
explained below. By inspection there is good point by 
point consistency between the present data and those 
of Ref. [i^ where the statistical accuracy was much 
higher, but where the temperatures studied did not quite 
span the critical temperature. The consistency between 
the two data sets implies that full equilibrium has been 
reached. The data in Figs. [T] and [2] for fixed (3 should 
tend to straight lines at criticality, curving upwards and 
downwards respectively at large L for P higher and lower 
than Pc', one can estimate Pc = 0.505(1) from this crite¬ 
rion. The value is mar gina lly higher than the estimate 
Pc = 0.5023(6) of Ref. [15j. Data we have obtained on 
other dimensionless parameters lead to very similar esti¬ 
mates for Pc 0. 

It is important to consider the value of the irrelevant 
scaling field thermal correction exponent 6 or equiva¬ 
lently the finite size correction exponent lo = Ojv. The 
first term in the RGT e-expansion for the ISG leading 
irrelevant operator is 0{d) = (6 — d) [ 130 , so 0(3) « 3, 
0(4) Ri 2, 0(5) ~ 1 (see [iJ] for the analogous site per¬ 
colation e-expansion). These values are only indicative 
as they are obtained from the leading terms of a series 
which (like the e-expansion series for the other exponents 
in ISGs) has never been summed. However, the values 
are qualitatively consistent with the observed 9{d) and 
uj{d) = 0{d)/h'{d) values for d = 3, 4 and 5. FSS esti¬ 
mates in 3d are lo = 1.12(10) and v = 2.56(4) [l3| so 
0(3) = uji^ ~ 3. HTSE estimates for various ISG mod¬ 


els in 4d are 0(4) Ri 1.5 and in 5d, 0(5) r; 1.0 [20L l21|. 


Simulation estimates in 4d and 5d are consistent with 
these values [l3113113 • Potentially there can also be an 
analytic correction with exponent Wa ~ > 2 in 

4d. If an analytic contribution is indeed present in 4d in 
addition to the leading conformal correction, the effec¬ 
tive exponent of the combined correction will be slightly 
increased. 

The relative strengths of the leading conformal cor¬ 
rection and the analytic correction will depend on the 
observable so the effective exponent for the combination 
can be expected to vary somewhat from observable to 
observable within a narrow range of values. 

The leading irrelevant scaling field correction is by def¬ 
inition the correction of this type which has the smallest 
exponent 0. From the overall agreement between the 
e-expansion and observed values, one can be fully con¬ 
fident that the leading irrelevant operator correction in 
4d is a;(4) r; 1.5, and hence no hypothetical correction 
term with a much smaller exponent (which if it existed 
would modify the behavior at very large L) can exist. 
Extrapolations of observed FSS data to infinite L using 
a correction exponent w Ri 1.3 are valid, which justifies 
the natural straight line extrapolations at Pc in Figs. [1] 
and m Similarly valid fits to the ThL data in Figs. [5] 



FIG. 3: (Color online) The 4d bimodal ISG. Thermodynamic 
derivative peak Umax = [9I/(/3,1/)/d/3]max data for dimen¬ 
sionless observables U = Wq,h,g (black squares, blue circles, 
red triangles). The peak location PmaxiL) against the inverse 
peak height l/Dmax(T). 


and [7] are made with a leading correction term having 
an exponent 0 r; 1.5. No hypothetical small exponent 
correction term producing ’’reentrant” behavior exists. 

For each dimensionless observable U{P,L), extrapola¬ 
tion to infinite L at Pc gives an estimate for a critical pa¬ 
rameter U{Pc,oo) characteristic of the universality class. 
For the bimodal ISG in dimension four these data show 
g{Pc,oo) = 0.525(5) and {p/L){Pc,oo) = 0.485(5). 

With the FSS rule Eq. (|4]) for a dimensionless observ¬ 
able U{P, L) the critical exponent v can be estimated 
in principle from [dU{P, L)/dP)]i3^ = To obtain 

these u estimates, the fits are multi-parameter as they 
involve simultaneous estimates for Pc, and the param¬ 
eters A, B and C as well as v. Nevertheless in practice 
reasonably precise values for v can be obtained. From 
analyses of data sets such as those shown in Figs. [T]and 
[5] on a number of different dimensionless observables, a 
global estimate for the bimodal model is = 1.13(1) [l6|. 


THERMODYNAMIC DERIVATIVE PEAK 
ANALYSIS 


Near criticality in a ferromagnet for dimensionless ob¬ 
servables U{P, L) the heights of the peaks Dcna^U, L) of 
the thermodynamic derivatives D{U,L) = dU{P,L)/dp 
scale for large L as 0,0 

[dU{P, L)/dp]cna. CX (l + + • • • ) (5) 





4 



T 



FIG. 4: (Color online) The 4d bimodal ISG. The Binder 
cumulant derivative peak height maximum [9p(/3, L)/cI/3]max 
normalized by L, against L (log-log scaling). 
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FIG. 5: (Color online) The 4d bimodal ISG. The effective 
susceptibility exponent 7 (r, L) with Pc = 0.505. Sizes L = 14 
to 5 from left to right. The red continuous curve is calcu¬ 
lated from the tabulated HTSE terms in Ref. [^. The green 
continuous curve is an overall fit. 

The temperature location of the derivative peak 
l3max{U,L) similarly scales as 

13c - /3max(C/, L) cx ^-^/"(l + + • ■ •) (6) 

As both /3c-/3max(t/, L) and l/T>max(C^, L) vary as 
at large L, a plot of the peak locations against the inverse 
peak heights tends linearly to Pc at large L. These esti¬ 
mates of Pc are independent of the FSS estimates. The 
observables used for U (/3, L) [l^ can be for instance the 


FIG. 6: (Color online) The 4d bimodal ISG. The effective 
correlation length exponent v{t, L) with Pc = 0.505. Sizes 
L = 14 to 4 from left to right. The green continuous curve 
is an overall fit. The red arrow is the exact bimodal high 
temperature limit i/(r = 1) = (d — 1/3)/?^. 



^(P,L)/P 


FIG. 7: (Color online) The 4d bimodal ISG. The ratio 
xiP,L)/{P.{P,L)/P)^ as function of ^{P,L)/P (log-log scal¬ 
ing). Sizes Z/ = 6, 8,10,12,14 from left to right. The green 
continuous curve is an overall fit. 


Binder cumulant g(P, L) or the logarithm of the finite 
size susceptibility ln(x(/?,L)). 

For an ISG just the same thermodynamic differential 
peak methodology can be used as in the ferromagnet. As 
far as we are aware this analysis has not been used pre¬ 
viously in the ISG context. In Fig. [3] we show data for 
peak height locations /3max(G, L) against the inverse peak 
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height 1/D^a.x{U,L) for the derivatives dU{j3,L)/dj3 
with dimensionless observables the Binder cu- 

mulant g{(3,L), together with Wq{(3,L) and h{l3,L) de¬ 
fined by 


and 


Wq{P,L) = 


h{P,L) = 


TT — 2 


m] 


- 2 


[(g2)]3/2 


-78 


(7) 


( 8 ) 


where the coefficients have been chosen so that the pa¬ 
rameters go from 0 at high temperature to 1 at low tem¬ 
perature. In Fig. 13] the linear extrapolations of the three 
sets of data points lead consistently to /3c = 0.505(1), in 
full agreement with the FSS data. 

From the thermodynamic derivatives, ’’The critical 
exponent v can be estimated without any considera¬ 
tion of the critical coupling /3c” 1^. As the exponent 
ly is close to 1 in 4d, we choose to plot data for the 
normalised Binder cumulant derivative ln(Di„ax(A)/L) 
against ln(L), Fig. [d] Corrections to scaling are visible 
for small L, but for L > 4 the slope is —0.0890(5), which 
corresponds to ty = 1.13(1), in full agreement with the 
FSS estimate above. 


THERMODYNAMIC LIMIT DERIVATIVES 

Analyses using the standard RGT scaling variable 
t = 1 — T/Tc = 1 — /3c//3, either in ferromagnets or in 
spin glasses, are restricted to the critical region and the 
FSS regime L ^ because t tends to diverge at 

high temperatures so scaling corrections automatically 
proliferate outside the critical regime. This is unneces¬ 
sary. A natural scaling variable, which has been used for 
the scaling of the susceptibility in ferromagnets for more 
than 50 years, is r = 1 — /3//3c The Wegner ther¬ 

modynamic limit (ThL) susceptibility scaling expression 
(^ . which is expressed in terms of r and is valid from 
criticality to infinite temperature,is 


x(t, oo) = C^T (I -h + h^T^ -) 


(9) 


At criticality r becomes identical to t and tends to 1 
at infinite temperature, so the corrections in the Weg¬ 
ner expression are well behaved for the entire paramag¬ 
netic temperature range. The first correction term in 
the equation is the leading confluent correction and the 
second an analytic correction. There is a closure condi¬ 
tion as x(r = 1) = I. In ISGs as the interaction energy 
is (J^), not (J), and the appropriate scaling variable is 
T = 1 — (/3//3c)^ or Tu, = 1 — (tanh(/3)/tanh(/3c))^. This 
is the scaling variable which was used from the beginning 
in ISG HTSE work 2^, 2^ [28143^ but curiously not in 


most analyses of numerical simulations. The Wegner ex¬ 
pression for the susceptibility [1^ carries over, with the 
appropriate definition for r. 

It has been pointed out that in Ising ferromagnets the 
ThL expression for the second moment correlation length 
^ (r) analogous to the Wegner expression for the suscep¬ 
tibility is 0, [ 23 ,0,131 1 


e(T) = Q/3i/v-‘' (1 + 


-I- ■ 


') 


( 10 ) 


The factor /3^/^ arises because the generic infinite tem¬ 
perature limit behavior is —>■ 1. 

In ISGs the factor /3^/^ becomes /3, again because (J^) 
replaces (J) [sT 


so 


ar) = (1 + 0^7 + • • •) 


( 11 ) 


with the ISG r. Following a well-established protocol 
27 , 3^ one can define a temperature dependent effective 
exponent for the susceptibility 

7(7 = -ain(x(T))/ain(T) (12) 

and for the second moment correlation length 

i/(t) =-51n(7T)//3)/ain(T) (13) 

which tend to the critical 7 and v respectively at criti¬ 
cality, and to and {d — K 1^)131 respectively in the 
high temperature r = 1 limit, where K is the kurtosis of 
the interaction distribution. {K = 1 in the bimodal case 
and K = 3 for the Gaussian). 

For data in the ThL regime, the Wegner expression 
with two correction terms Eq. m translates exactly into 


7(r) = 7 - 


1 -I- + b^ry 


(14) 


where the second term is an effective subleading correc¬ 
tion. There is an equivalent relation for iy{T). Fixing 
/3c = 0.505, the 7 (r) and plots for all the different 
sizes L are shown in Fig. [5] and Fig. | 6 l In Fig. [5] the 
effective 7 (/ 3 ) evaluated from an explicit sum of the first 
15 terms in the HTSE series of Ref. 2l| is also shown. 

The fit curves are adjusted to those data which are 
in the ThL regime, and the fit is then extrapolated to 
criticality to obtain estimates for the critical 7 and ly. 
Just as for uj in the FSS regime, 9 is the lowest correction 
exponent and so a further hypothetical correction term 
with a much lower exponent which could modify the form 
of the fit curve between the ThL data and criticality can 
be ruled out. The 7 (t) and ;z(t) plots are however rather 
sensitive to the exact value of /3c. 

The fits shown correspond to an ISG temperature de¬ 
pendent susceptibility effectively fitted by 

x(/3) = 0.634r-2 ™ (l -k 0.38 t^-^ -k 0.227 - 0.022 t®) 

(15) 
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and a second moment correlation length 


C(/3) = 0.914/3r"^-^^ (l + 0.055 t^-^ + 0.042r3) (16) 


so with estimated critical exponents 7 = 2.70, v = 1.12, 
Q = 1.5, with a leading correction term and further higher 
order terms. 

The ratio 


v{t,L) 


^lnx(/?^T)/^lnr 
d ln(^ (/ 32 , L) //3)/9 In r 


2-7y(/3,i) (17) 


does not involve /3c, so r]{l3,L) can be measured with¬ 
out any knowledge of /3c. Then the ri{j3, L) as func¬ 
tions of /3/^(/3, T) can be extrapolated to /3/^(/3,L) = 0 
for infinite L to estimate the critical value of 77 with¬ 
out involving /3c. Alternatively we plot, see Fig. [3 
y = x(/3,i)/(C(/?,i)//3)^ against x = \{P,L)/j3. At in¬ 
finite temperature y = x = 1, and towards criticality 
for the ThL regime data y oc and 

X oc T~'^ . Hence the slope of ln( 2 /) against ln(a:) tends to 
—rj as criticality is approached. The fit and limiting slope 
can be estimated without specifying an explicit value of 
Pc- A satisfactory fit curve with a single effective correc¬ 
tion term is y = Cx~'^{l {C~^ — l)x~^) with rj, 9 and 

C adjustable parameters. The values in the fit shown are 

= -0.415(10), 9 = 1.3, and C = 0.78. 

The critical exponents estimated without any knowl¬ 
edge of Pc are therefore v = 1.13(1) from the thermody¬ 
namic derivative result Fig.0]and 77 = —0.415(10), and so 
7 = 17(2 — 77 ) = 2.73(3) combining these values. The val¬ 
ues are in excellent agreement with the estimates from 
the extrapolated ThL derivatives of Fig. [S] and Fig. [S] 
which validates the Pc estimate and the overall method¬ 
ology. 


HIGH TEMPERATURE SERIES EXPANSION 



1 /{n+1) 


FIG. 8: (Color online) The 4d bimodal ISG. The HTSE ratio 
series R{n -|- 1) for the higher order susceptibilities r 4 (red 
squares) and Fa (blue circles). Green curves : fits. Arrow : 
critical point. 


with T = 1 — w/wc, the corresponding HTSE terms are 


i{n) = 


Cn'^- 


1 -h 


r(7) 


T{'~f)wc V ' r(7 — 0) n® 
The ratios between successive terms are 
a{n -|- 1 ) 


R{n + 1) = 


a{n) 

7-1 


= — ( 1 + 

Wc \ n + l 


1 - 


r(7) 


9a^ 


r(7- 6») (n-h 1)®+1 


( 20 ) 


( 21 ) 


In the HTSE approach, series of exact terms are eval¬ 
uated, which are summed to obtain the temperature de¬ 
pendence of the susceptibility. In ISGs the series for the 
spin glass susceptibility is 


x{w) = 1 -b 0 ( 1 ) 77 ; -b a(2)w^ + a(3)w^ -b ■ 


where w = tanh^(/3 


(18) 


201 . |29j , or the same equation with 
w replaced by /3^ [2l|. The terms a(n) are exact. The 
number of terms which can be calculated is limited by 
practical considerations, and up to now tt, = 15 has been 
an upper limit in ISGs. The series have been analysed 
by Dlog-Pade, Ml and M2 techniques [iSllllii which 
are not particularly transparent for a non-specialist. An 
alternative classical approach is the ratio method [3, [ 23 , 
33|. For a model where the temperature dependeuce of 
the thermodyuamic limit (infinite size) susceptibility is 

in 


In favorable cases the ratio series R(n + 1) plotted as 
a function of l/(n -b 1 ) can be extrapolated to infinite 
n to estimate the critical temperature from the inter¬ 
cept, 1/wc, the critical exponent 7 from the initial slope 
dR(n + l)/d(l/(n + 1)) = (7 — l)/77;c, and the confor¬ 
mal correction exponent 9 from the leading non-linear 
term. Unfortunately there can be parasitic oscillating 
terms in the ratios, arising from anti-ferromagnetic poles 


27| . particularly in simple hypercubic lattices. It turns 


X(w) = C^T ^ (1 -ba^r® + ••■) 


(19) 


out that in the ISGs the parasitic terms in the suscep¬ 
tibility series become rather strong when the dimension 
drops. Explicit lists of terms for r 2 in dimension <3 = 4 
and below were tabulated in Ref. [31 • 

For the simple hypercubic bimodal ISGs iu general di¬ 
mension Klein et al. (3| tabulated the terms to order 
15 not only for the ISG susceptibility x(^) series (called 
r 2 in their nomenclature) but also for the higher order 
susceptibilities F 3 and r 4 which they dehne. These sus¬ 
ceptibilities have critical exponents 73 = {3x + diy)/2 and 
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74 = 27 + dv respectively [20|, [27|. Unfortunately the 


parasitic terms are so strong in the 4d r 2 susceptibility 
series that the ratios cannot be readily exploited so as to 
make accurate estimates of the critical parameters. This 
explains why the error bars quoted in Ref. 2l| for 
and 7 in 4d are rather large. However, we have summed 
the tabulations of for Fa and r 4 to obtain the 03 ( 71 ) 
and 04 ( 71 ) for these series, and have found that for these 
higher order susceptibilities the ratios behave very regu¬ 
larly in dimension 4 (and above), as can be seen in Fig.| 8 l 
Expression EH can be applied to these ratios and reliable 
extrapolations can be made through which one can es¬ 
timate TCc, 7 ( 3 ), 7 ( 4 ) and 6. Optimal fits excluding the 
small n values are 


, , f 5-3 4.0 \ 

i? n-Pl =4.60 1-f - - --7 - — 5 ^ 22 

^ ' \ (n + l) ( 77 -fl) 2 - 35 y ^ ' 

for Fa, and 

Rin + 1) = 4.60 (^1 -f (23) 


for F 4 . The joint intercept 1/wc = l/tanh^(/3c) = 
4.60(2) corresponds to /3c = 0.505(2), and the initial 
slopes correspond to 73 = (87 -I- di^)/2 = 6.30(10) and 
ji = 2^ + dv = 10.15(10). The leading correction terms 
have an exponent 9 « 1.4, and prefactors with nega¬ 
tive sign. By inspection any hypothesis of another term 
(which would then be the true conformal correction to 
scaling) having a much lower exponent and with a pref¬ 
actor of the opposite sign can once again be ruled out. 
There are no statistical errors as all points are exact; sys¬ 
tematic errors arising from the extrapolation are small. 

The exponents estimated above from the simulations 
and so quite independently from the HTSE data corre¬ 
spond to 73 = (37 -I- 4i/)/2 = 6.3(1) and 74 = 27 -|- 4i7 = 
9.92(15), and so are fully consistent with the HTSE val¬ 
ues. 



1 / 1 '-" 


FIG. 9: (Color online) The 4d Gaussian ISG. FSS for the cor¬ 
relation length ratio ^{j3,L)/L at fixed inverse temperatures 
P = 0.5625, 0.5600, 0.5575, 0.5550, 0.5525, 0.5500 from top 
to bottom, against 1/. Dashed line : critical behavior. 


ca 
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GAUSSIAN INTERACTION MODEL 


We will present equivalent data for the hypercubic ISG 
model with Gaussian interactions in dimension four. All 
the discussions above apply equally well to the Gaussian 
data. Corrections for the Gaussian model are in general 
much weaker than for the bimodal model. Unfortunately 
no HTSE higher order susceptibility data are available 
for the Gaussian model. 

The 4d Gaussian critical inverse temperature from 
both the FSS data. Figs. [9] and [TOl and the thermo¬ 
dynamic derivative peak locations Fig. [TTl is estimated 
to be Pc = 0.5555(10) or Tc = 1.800(3). This value is 
fully consistent with previous FSS estimates from Binder 
parameter and correlation length ratio measurements 

0.555(3) [11,133 and 0.554(3)0. 


FIG. 10: (Color online) The 4d Gaussian ISG. FSS for the 
Binder cumulant g{P, L) at fixed inverse temperatures /3 = 
0.5625, 0.5600, 0.5575, 0.5550, 0.5525, 0.5500 from top to 
bottom, against \/Ld'^. Dashed line : critical behavior. 


The Gaussian critical parameters for the normalised 
second moment correlation length p{P,L)/L and the 
Binder cumulant (?(/3, L) from Figs. [9landfT0l are {^/L)c = 
0.458(5) and gc = 0.485(5) respectively. The value for the 
Binder cumulant is close to the estimate gc = 0.470(5) 
given in Ref. 0 . 

The exponent v estimated directly form the Binder cu¬ 
mulant thermodynamic derivative peak heights, Fig. 1121 
is 1.030(5). This value is fully consistent with the esti- 
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FIG. 11: (Color online) The 4d Gaussian ISG. Thermody¬ 
namic derivative peak Dmax = L)/9/l]max data for 

dimensionless observables U = Wq,h,g (black squares, blue 
circles, red triangles). The peak location /3max(T) against the 
inverse peak height l/Ilmax(T). 



FIG. 12: (Golor online) The 4d Gaussian ISG. The Binder 
cumulant derivative peak height maximum [clgi((l, L)/cl/3]max 
normalized by L, against L (log-log scaling). 

mate iz = 1.02(2) given in Ref. 

The plot of ln(?/) against ln(a;) with y = 
xil3,L)/kil3,L)/l3)‘^ andx = ^{I3,L)/I3, Fig. [131 is fitted 
with a single effective correction : 

y{x) = 1.2a;°-305 (l _ 0.17x“^-^) (24) 

corresponding to a critical exponent y = —0.305(10). 
This value is again fully consistent with but more ac¬ 
curate than the estimate rj = —0.275(25) given by 


FIG. 13: (Golor online) The 4d Gaussian ISG. The ratio 
X{P,L)/{£,{P,L)/Pf as function of ^(d,i)//3 (log-log scal¬ 
ing). Sizes L — 6,8, 12 from left to right. 



T 


FIG. 14: (Golor online) The 4d Gaussian ISG. The effective 
exponent 7 (r) with /3c = 0.5555. Sizes L = 12,10, 9, 8, 7, 6, 5 
from left to right. The blue continuous curve is calculated 
from the tabulated HTSE terms in Ref. [^ . The red arrow 
is the exact high temperature limit 7 (r = 1) = 2ddc- 

The fit to the ThL effective 7 (r) with a fixed /3c = 
0.5555, Fig. [T31 corresponds to 

X = l.Olr-2-3^ (1 - 0.010ri°) (25) 

so with a critical exponent 7 = 2.37(2) and only a very 
weak high order correction term. The prefactor for the 
usual leading correction term appears to be accidentally 
very close to zero. A similar conclusion was reached from 
the FSS analysis of the 4d Gaussian model in The fit 
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FIG. 15: (Color online) The 4d Gaussian ISG. The effective 
exponent v{t) with /3c = 0.5555. Sizes L = 12,10, 8, 7, 6, 5, 4 
from left to right. The red arrow is the exact Gaussian high 
temperature limit = 1) = (d — l)/3c. 


to the ThL effective i^{t) with fixed Pc = 0.5555, Fig. [T51 
corresponds to 

C(/3) = (1 + O.OTOr^ ®) (26) 

so a critical exponent i/ = 1.025(3) and a correction term 
with a similar exponent to the exponent of the bimodal 
model leading correction. 

The critical constants and exponents for a 4d diluted 
bimodal model Ref. @ were estimated to be g{Pc,oo) = 
0.472(2), (^/L)(/3c,oo) Ri 0.440, 7 = 2.33(6), and jz = 
1.025(15). Each of the diluted bimodal critical values is 
similar to that for the 4d Gaussian model, but the diluted 
bimodal model and the standard bimodal model studied 
above have quite different critical properties. 


CONCLUSION 

We compare in TableUthe bimodal and Gaussian inter¬ 
action distribution model estimates for each of the critical 
parameters and critical exponents studied above. For all 
the parameters and exponents the estimates for the two 
models are quite different, with in each case differences 
of the order of 10%. We conclude that in 4d ISGs the 
critical exponents, like the values of the critical ordering 
parameters, depend on the form of the interaction dis¬ 
tribution. ISGs in dimension four with different interac¬ 
tion distributions are clearly not in the same universality 
class. 


TABLE I: Values of the 4d bimodal and Gaussian critical 
parameters and exponents. 


Parameter 

Bimodal 

Gaussian 

Pc 

0.505(1) 

0.5555(10) 

giPc,oo) 

0.525(5) 

0.485(5) 

(C/L)(/3c,oo) 

0.485(5) 

0.458(5) 

7 

2.70(3) 

2.37(2) 

V 

1.113(1) 

1.030(5) 

V 

0.42(2) 

0.305(10) 
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